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Based on a sample of 58 x W 6 J/ip decays collected with the BESII detector at the BEPC, the A 
decay parameter for A — ¥ piv + is measured using about 9000 J/if) — > AA — > ppTv + ir~ decays. A 
fit to the joint angular distributions yields a^(A —¥ piv + ) = —0.755 ± 0.083 ± 0.063, where the first 
error is statistical, and the second systematic. 
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I. INTRODUCTION 

Nonleptonic hyperon decays were very important in early studies of parity violation in particle physics [l| ■ In 
a hyperon decay, Y —> Bit (y:hyperon, B:baryon), the angular distribution of the baryon in the center-of-mass 
(CM) system of Y takes the form 4^ oc 1 + ayPy ■ Pd, where Py is the polarization vector of the hyperon, pd is 
the momentum unit vector of the baryon and ay is the hyperon decay parameter, which characterizes the parity 
violation in hyperon nonleptonic decays. For the decay A —> pir~ , the decay parameter is a a = 0.642 ±0.013 [2]; 
for the counterpart A — > pir + , the decay parameter was only measured by the DM2 Collaboration with large 
uncertainty: a A = —0.63 ±0.13 p|. 

A precise measurement of the A decay parameter allows a more sensitive search for CP symmetry violation 
in A decays 0, A CP-odd observable, A, may be defined as 

A= aA + a \ (1) 
ola - a A 

If CP is conserved, a\ = —olr an d this observable vanishes; while a nonzero value of A implies evidence for 
CP asymmetry in A decays. Previous analyses to search for CP asymmetry in A nonleptonic decays have been 
performed at pp colliders by the R608 @ and PS185 0] Collaborations, and at an e + e~ collider by the DM2 
Collaboration [3j, but the precisions of the measurements are limited by low statistics. 

A precise A decay parameter is also essential in determining the Q + or S + decay parameters. Nonpolarized 
fi + or S + decays can produce polarized A particles, so in the A rest frame, the angular distribution of the final 
state antiproton, takes the form d ^ g oc 1 + a^a\ cos 9 To extract from the product a^a^ the value of 
a A is required. A similar argument holds for S + decays. 

Pair production of AA in J/tp decays allows one to study the A decay parameter. Although the A and A 
particles are nonpolarized in J/ip decays, their helicities are correlated by helicity conservation. Hence, the A 
decay parameter can be extracted from the helicity correlation between A and A. Experimentally, this decay 
is very clean and can be reconstructed with high efficiency by selecting events with four charged tracks. The 
58 x 10 J/tp decays used for this analysis were taken with the BESII detector at the BEPC storage ring at a 
center-of-mass energy corresponding to Mji^. They offer an opportunity to measure a more precise A decay 
parameter. 



II. THE BES EXPERIMENT 



The Beijing Spectrometer (BES) detector is a conventional solenoidal magnet detector that is described in 
detail in Ref. [9(; BESII is the upgraded version of the BES detector [Toj]. A 12-layer vertex chamber (VTC) 
surrounding the beam pipe provides trigger and track information. A 40-layer main drift chamber (MDC), 
located radially outside the VTC, provides trajectory and energy loss (dE/dx) information for charged tracks 
over 85% of the total solid angle. The momentum resolution is cr p /p — 0.017\/l + p 2 (p in GeV/c), and the 
dE/dx resolution for hadron tracks is ~ 8 %. An array of 48 scintillation counters surrounding the MDC 
measures the time-of-flight (TOF) of charged tracks with a resolution of ~ 200 ps for hadrons. Radially outside 
the TOF system is a 12 r.L, lead-gas barrel shower counter (BSC). This measures the energies of electrons and 
photons over ~ 80% of the total solid angle with an energy resolution of cte/E = 22%/y/E (E in GcV). Outside 
of the solenoidal coil, which provides a 0.4 T magnetic field over the tracking volume, is an iron flux return that 
is instrumented with three double layers of counters that identify muons of momentum greater than 0.5 GeV/c. 

A GEANT3 based Monte Carlo (MC) program with detailed consideration of detector performance (such as 
dead electronic channels) is used to simulate the BESII detector. The consistency between data and Monte Carlo 
has been carefully checked in many high purity physics channels, and the agreement is quite reasonable [Tlj ]. 
Several J/tp — > AA MC samples are generated and used for determining the detection efficiency, determining 
the normalization factor in background subtraction, and performing an input-output check. 



III. EVENT SELECTION 



The A is reconstructed under the assumption of A — > pix~ decay, and A under the A — > pir + hypothesis. A 
candidate track is required to have a good helix fit, the polar angle must satisfy | cos# c h| < 0.8, and transverse 
momentum p xy > 0.07 GeV/c. Events are required to have four charged tracks with total charge of zero. Protons 
and antiprotons are identified using TOF and dE/dx information with the requirement that their confidence 
levels be larger than 0.01. A four constraint (4C) kinematic fit is applied under the pftn + 'K~ hypothesis, and 
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FIG. 1: Distributions for selected J/ip —¥ AA candidates, (a) xlc- Dots 
with error bars are data; histogram is for MC simulation, (b) Scatter 
plot of m p7r - versus m pv +. 



xlc < ^0 * s required. A comparison of xlc distributions of data and MC simulation is shown in Fig. Q] (a). 
The scatter plot of mp n + versus m p7r - is shown in Fig. [T] (b), and J/ip — >• AA is clearly seen. The A and A 

are selected by requiring |M p7r - - 1.1156| < 0.015 GeV/c 2 and \Mp„+ - 1.1156| < 0.015 GeV/c 2 . The mass 
distributions of pir^ and pTr + are shown in Fig. [5J 
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FIG. 2: Comparisons of invariant mass distributions between data and MC. Histograms are 
normalized MC; points with error bars are data, (a) m p7r - and (b) rrip w +. 



After applying the above selection criteria, 8997 J/ip — > AA events are selected. The ang ular distribution of 
A in the J/ip rest frame is found to be consistent with the dedicated analysis in Ref. [12] within uncommon 




statistical errors. 

Backgrounds are studied with MC simulations. The main backgrounds come from J/tp — > E°S° — > 2jAA — > 
2jpir~pTT + (28 events), J/ip AE° — > 7AA — > ~/pTr~pTr + (36 events), J/ip — > AT, — > 7AA -> jpTr~pir + 
(37 events), <//?/> — > p7r~p7r + (10 events), J/V 1 — > A ++ A (5 events), J/V-> — > A ++ p7r~ (1 events), and 
J/ip —> A p7r + (8 events). The fraction of remaining background events in the data sample is about 1.4%. 
Background decays to intermediate states including A or A are generated with the helicity amplitude method 
[lj|, in which hyperon decay parameters are set according to PDG values Q. The effect of the background 
contamination on the physics results is included as one source of systematic error. 



IV. FIT TO DATA 



The angles used in this analysis are defined in the helicity frame, shown in Fig. [3] For J/tp — > AA, the z axis 
of the J/tp rest frame is along the A out-going direction, and the solid angle flo(9, 4>) is between the e + direction 
and the A out-going direction. For A — >• pir~ , the solid angle of the daughter particle f2,-(0j,0j) is referred to 
the A rest frame, and the z axis is also the A out-going direction. For A we use symbols with bars. 

For J /if) — > AA — > pn~pn + , the partial decay rate is [13[ 

oc (1 — a) sin 2 8 [I + a\a^(cos9i cos#i + sin6>i sin#i cos(0i + (pi))] 

- (1 + a)(l + cos 2 6)(a\ax cos^i cos^i - 1), (2) 

where dfl — dflgdClidQi, a is the angular distribution parameter for A, and oa^a) is the A (A) decay parameter. 
Equation. Q only allows one to obtain the product oa^a from a fit to data. To extract the value of a a, one 
needs to fix the A decay parameter a\, say, at the world average value cha = 0.642 2]. To determine the 
CP-odd observable defined in Eq. {TJ, one may make a replacement 



"A^A = ]T^7 a A ( 3 ) 



in Eq. ©. 

An unbinned maximum likelihood method is used to fit the data. As widely used in partial wave analyses (l4|. 
a normalized probability density function for J/ip — > AA — > p-K~pn + is defined as 

Prob(f!,a,a A ) = ^^, (4) 
a 

where f2 = (f2o,^i,^i), and a and are parameters to be determined. da/dQ is given by Eq. a is the 
total cross section given by 



(5) 
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FIG. 4: Comparison between data (dots with error bars) and fit results (histograms), (a) Distribution of 
cos6> of A in J/ip rest frame, (b) Distribution of cos(9i cos#i + sin^i sin#i cos(<j!>i + 4>\). 



where e(f2) is the detection efficiency. The total cross section a can be determined by MC numerical integration 
with a phase space generator over the allowed kinematic region of J/rp — > AA — > pTt~piT + , i.e. 



= 1 y C ( da 
a ~ /V, 1 WO 



N M c frf [dfl 

where Nmg is the number of selected MC events. 

The maximum likelihood function is given by the joint probability density of the selected J/t/j — > AA 
pTT~pTT + events 



(6) 



N Bvt 

C = II Prob(r2,a,a A ), 
i=i 



(J) 



where N evt is the number of selected data events. To determine the unknown parameters a and a^, the function 
S, which is defined as 



S = -ln£, 

is minimized. 

To check the goodness of fit in our analysis, we define 



(8) 
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(9) 



over the scatter plot of cos0 versus | cos#i cos$i ±8^6*1 sin(9i cos(</>i + Here N is the number of cells, nf r 
and v!f c are the numbers of events in the ith cell of the scatter plot for data and MC simulation, respectively. 
Such a variable should be distributed according to the % 2 distribution with ndf — N — K degrees of freedom, 
where K = 2 is the number of parameters to be determined in our fit. In this experiment, data and MC events 
are divided into 8 x 10 cells in the scatter plots. 
Fitting the data using Eq. ([8]), we obtain 



0.70 ±0.06, 



(10) 



a K = -0.755 ±0.083, 



(11) 



G 



while fitting the data using Eq. ([SJ with parameters a and A, we obtain 

A = -0.081 ± 0.055, (12) 

where the errors are statistical only. Comparisons between the data and the fit results are shown in Fig. 2] (a) 
and (b). The fit yields x 2 /ndf = 78.59/(80-2) = 1.01. If we force A = 0, we obtain a A = -a A = 0.696±0.038, 
where the error is statistical only. 

V. INPUT-OUTPUT CHECK 

To validate the fitting procedure, a MC sample of 2 x 10 6 J/tp — > AA — » ppir + ir~ events is produced according 
to Eq. @. The input parameters are a = 0.62 and a\ = —a A = 0.642. The MC sample is required to pass 
the same selection criteria as used for data selection, and the same fitting procedure is applied to the selected 
events with a A fixed to 0.642. The fit yields a = 0.612 ± 0.010 and a A = -0.640 ± 0.013, consistent with the 
input values to within the la statistical errors. 

VI. SYSTEMATIC ERRORS 

A. Background contamination 

Contamination from background channels is studied using MC samples, including J/tp — > £ S°, AE° + 
c.c., A ++ A , A ++ p-K~, A~~p7T + and ppn + n~ . Background channels with A/A intermediate states are 
generated with full helicity amplitude information, and the decay parameters for A, S° and their antiparticlcs 
are fixed at the world average values |2[ . These samples are required to pass the same selection criteria as used 
for the the signal channel J/tp —> A A —> ppir + ir~~ . The selected background events are subtracted from the data 
in the fit. The difference of the A decay parameter a A from the fit without background events, 5a A = 0.021, is 
taken as the uncertainty from background contamination. 

B. Decay parameter a A 

The systematic error due to the uncertainty on the A decay parameter [2j is determined. In fitting to data, 
the central value of a A is obtained by fixing a A = 0.642; the difference from the central value, 5a A = 0.015, is 
obtained by changing a a by 1 standard deviation (±0.013). 

C. MC simulation and detector response 

The consistency between data and the MC simulation of the detector response for J/tp — > AA events can 
be determined using the channel with the same final states J/tp — > ppir + ir~ after rejecting A, A, A ++ and 
A~~ intermediate states. Therefore fitting to this sample with Eq. © should yield aACt\ — 0. The difference 
from zero gives the systematic error due to the MC simulation of the detector response. In fitting this sample 
with Eq. angles for A or A are replaced with ones for the quasi-two body systems (ptt~ and pir + ). The 
fit yields aACt\ — 0.008 ± 0.036; the uncertainty from MC simulation and detector response is taken as 0.044, 
the linear sum of the deviation and its uncertainty. The resulting distributions of cos 9 and cos 9\ cos 9\ + 
sin 9\ sin 0i cos(</>i + ipi) are shown in Fig. [5] 

D. Hadron interaction model and wire resolution 

The systematic error associated with MC simulation of hadronic interactions of final state particles with the 
detector is estimated with two models, GCALOR and FLUKA [l6j]. The central value is obtained with 
GCALOR model, while the difference from the central value using FLUKA model, 8a A — 0.005, is taken as the 
systematic error. The uncertainty associated with the simulation of the MDC wire resolution is estimated with 
two different versions of MC simulation software, which yields 6a A = 0.037. 

Table U summarizes all systematic errors. The total systematic error is estimated to be Sa A — 0.063, assuming 
all the sources are independent and summing them in quadrature. 
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FIG. 5: Fitting results of the control sample J/ip —> ppir + ir~ used for determining the uncertainty associated 
with the MC simulation and detector response. Points with error bars are data, and the histograms are the 
fit results, (a) Distribution of cos (9. (b) Distribution of cosSi cos Si + sin Si sin Si cos(<^i + <j>i). 



TABLE I: Systematic errors in a A . 



Backgrounds 0.021 
a A 0.015 
MC simulation and detector response 0.044 
Hadron model 0.005 
Wire resolution 0.037 
Total 0.063 



VII. RESULT AND DISCUSSION 

Including systematic errors, we obtain 

a A = -0.755 ± 0.083 ± 0.063, 
A= -0.081 ± 0.055 ± 0.059, 

where the first errors are statistical and the second systematic. The comparison between DM2 Q and our 
results for a and A is given in Table [TXJ They agree with each other within errors, while BES has improved 
precision. 

Our measurement gives the value of CP-odd variable A = —0.081 ± 0.081. The precision of this measurement 
is insufficient to observe CP violation at the level predicted by the standard model: A = —2.10 x 10~ 5 in the 
Kobayashi-Maskawa model or A = —1.10 x 10~ 4 in the Weinberg model [3(. The precision of this measurement 
is expected to be highly improved at BESIII in the near future with 10 10 J ftp decays accumulated. 
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TABLE II: Comparison of a^(A — > p-7r + ) and A with DM2 measurements 0]. 





«a(A ->ptt + ) 


A 


DM2 


-0.63 ±0.13 


0.01 ±0.10 


This work 


-0.755 ± 0.083 ± 0.063 


-0.081 ±0.055 ±0.059 
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